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The ultrafast dynamics of surface electromagnetic waves photogenerated on aluminum film per- 
forated with subwavelength holes array was studied in the visible spectral range by the technique of 
transient photomodulation with ~ 100 fs time resolution. We observed a pronounced blueshift of the 
resonant transmission band that reveals the important role of plasma attenuation in the optical re- 
sponse of nanohole arrays. The blueshift is inconsistent with plasmonic mechanism of extraordinary 
transmission and points to the crucial role of interference in the formation of transmission bands. 
The transient photomodulation spectra were successfully modeled within the Boltzmann equation 
approach for the electron-phonon relaxation dynamics, involving non-equilibrium hot electrons and 
quasi-equilibrium phonons. 
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Light transmission through periodic subwavelength ar- 
rays of holes in optically thick metal films JL,] (plasmonic 
lattices) has attracted intensive interest due to poten- 
tial applications in near-field microscopy, optoelectron- 
ics, and biosensing 0, On a smooth metal-dielectric 
interface light cannot couple to surface electromagnetic 
waves (SEW) such as, e.g., surface plasmon polaritons 
(SPP) because energy and momentum conservation can- 
not be simultaneously obeyed. But the periodicity in 
nanohole arrays leads to the formation of band struc- 
ture that enables light to directly excite SEW, which 
consequently mediate the optical transmission through 
the metal film. The transmission spectrum, To{u!), is 
then characterized by several resonant bands, each fol- 
lowed by a minimum caused by Wood's anomaly that is 
due to the diffraction grating related to the corrugated 
film [l| . These extraordinary optical transmission (EOT) 
resonances were originally attributed to grating coupling 
of light to SPP excitations on both sides of the perfo- 
rated film [U, 0, IB ! but later other possible explanations 
were also suggested such as interference of resonant cav- 
ity modes [a, Q and dynamical diffraction Qj, Q. The 
EOT spectrum has been shown to depend sensitively on 
geometric parameters such as film thickness, lattice pe- 
riodicity a, and holes size and shape as well as on 
the material optical constants such as substrate (sd) and 
hole (e/i) dielectric constants, and that of the metal film: 
em{<jj) = Em + *Em- Howcvei, relatively little is known 
about the EOT dynamics that arises when the optical 
constants are subjected to a pulsed time- dependent per- 
turbation. It has long been known 11| that time-resolved 
ultrafast spectroscopy can provide valuable insight into 
the makeup of underlying excitations from the transient 
changes of optical spectra [13 • 

Here we investigated the ultrafast dynamics of the 
EOT spectrum in Al-based plasmonic lattice measured 
by the transient photomodulation (PM) spectroscopy in 
the visible range with ~ 100 fs resolution. Surprisingly, 



we found a transient blueshift of the main EOT band, 
which unravels the important role of the plasma attenu- 
ation in determining the dynamical response of nanohole 
arrays. The blueshift is inconsistent with the simple 
SPP mechanism of EOT, and indicates the dominant 
role of interference in the formation of the transmis- 
sion bands. The obtained transient dynamics is mod- 
eled within the Boltzmann equation for electron-phonon 
relaxation with non-equilibrium hot electrons and quasi- 
equilibrium phonons, and subsequent heat transfer to the 
glass substrate. 

The plasmonic lattice sample was a 70 nm thick Al 
film, 5x5 mm^ in area deposited on a glass substrate. 
The metal film was perforated with circular holes having 
diameter D = 150 nm in a square hole array with lattice 
constant a = 300 nm resulting in a fractional aperture 
area of 25%. The ultrafast laser system used for measur- 
ing the transient PM spectrum was a Ti:sapphire regen- 
erative amplifier with pulses of 100 fs duration at photon 
energies of 1.55 eV, having 400 /xJ energy per pulse at 
a repetition rate of 1 kHz. The second harmonic of the 
fundamental pulses at 390 nm (3.1 eV) were used as the 
pump beam. The probe beam was a white light super- 
continuum generated in a 1-mm thick sapphire plate that 
covers the spectral range from 1.6 to 2.8 eV; the white 
light super-continuum was also used for measuring the 
unperturbed transmission, Tq(lu). The pump and probe 
beams were directed onto the perforated Al film from the 
air side to maximize coupling of light to the SEW. The 
pump and probe beams were directed onto the perforated 
Al film from the air side to maximize coupling of light 
to the SEW. The transient PM spectrum was obtained 
from the photoinduced change (AT) in Tq using a phase- 
sensitive technique with a resolution AT/Tq ^ 10"^ that 
corresponds to a photoexcited electron density in the Al 
metal of ~ 10^^ cm^'^ /pulse. 

Figure [1] shows the EOT spectrum, Tq{u}), of the Al 
plasmonic lattice at normal incidence. Several EOT 
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FIG. 1: (Color online) The measured (solid line) and calcu- 
lated (dashed line) normal incidence optical transmission, To 
through an Al plasmonic lattice comprised of a square array 
of nanoholes with lattice constant 300 nm and hole diameter 
150 nm. 



bands and transmission anti-resonances (AR) are ob- 
served. We note that To{u!) at resonance frequencies 
substantially surpasses the fractional aperture area of 
25%. The lowest frequency EOT band peaks at ~ 2.1 
eV; whereas the lowest frequency AR feature that is 
associated with Wood's anomaly at the Al/glass inter- 
face with wavevector fc^^'^) = g^gi^"^ (where Sgi = 2.43 
and g — Itx j a ) is in the form of a dip at 2.6 eV. 
The AR's at higher frequencies are higher-order Wood's 
anomalies that correspond to wavevectors = 

g{rr? +n^)^^^e^^^^ (m, n are integers) as follows: 3.7 
eV (fc(") = gSgl^'^i^l'^) that corresponds to the (1,1) 
branch at the Al/glass interface; and 4.0 eV (fc(°^^ = g) 
that corresponds to the (0,1) branch at the Al/air inter- 
face. It is noteworthy that the lowest frequency EOT 
resonance (at ~ 2.1 eV) occurs at considerably lower (by 
~ 0.4 eV) frequency than that predicted by the lowest- 
order SPP having wavevector that satisfies the Bragg 

condition Q: fci^/^ = g [(m^ -f n^) {e'-^ + e~^)\^'\ 

Figure [5] shows the PM transient response (photo- 
bleaching) at 650 nm (1.9 eV). Two distinct dynamic 
responses are apparent: (i) a PM rise during ~ 1 ps fol- 
lowed by a plateau [Fig.[2fa)] ; (ii) after about 3 ps the sig- 
nal starts recovering with characteristic time constant of 
~ 40 ps [Fig.[2Ib)]. These are unique features associated 
with the Al plasmonic lattice. In contrast to noble metals 
[ist , the interband transitions in Al do not play a signifi- 
cant role in determining e„j(cj), and therefore the PM re- 
sponse dynamics stems from the transient change in the 
electron-phonon (e-p) relaxation rate, jep(TL), where Tl 
is the metal lattice temperature. The dielectric function 
em(w) of Al is usually approximated by a Drude form; 
Em = 1 — ll>p/uj{lu + ij), where ujp is the bulk plasma 
frequency, and 7 = 70 + 7ep is the attenuation with rel- 
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FIG. 2: (Color online) Short-time (a) and long-time (b) dy- 
namics of the transient PM of the Al plasmonic lattice at 1.9 
eV. The solid line is a fit based on Eqs. ((2]|4]) in the text. 
Insets: (a) short-time evolution of the Al efi'ective electron 
temperature, T* (solid line) and lattice, Tl (dashed line); (b) 
long-time evolution of Tl ■ 



atively small contribution, tq, from electron-electron (e- 
e) and impurity scattering [14| . An increase in that 
is caused by energy exchange between the photoexcited 
(hot) electrons and phonons, is accompanied by an in- 
crease in 7ep(lL) due to larger phase space available for 
the phonon states at elevated temperature. The tran- 
sient change in em{i^) leads to the shift of EOT resonant 
frequencies (but less so for the AR frequencies, which are 
mainly determined by the plasmonic lattice periodicity). 
We thus interpret the PM rise of 1 ps as due to the 
increase in 7 (and consequently also in e^) caused by 
e-p energy exchange. Furthermore, due to the strong e-p 
interaction in Al, the excess energy transfer to the lattice 
occurs before the electron thermalization process is com- 
pleted; therefore the short-time dynamics is due to an 
inherently non-equilibrium process. After an equilibrium 
between the electron and phonon subsystems in the Al 
film is reached at higher Tl [plateau in Fig. [D^a)] , the ex- 
cess heat dissipates into the glass substrate within 40 
ps [Fig.Hb)]. 



3 



10 ps 
50 ps 
100 ps 



blueshift 




2.0 2.2 

Photon energy (eV) 



1.8 2.0 2.2 

Photon energy (eV) 



2.4 



FIG. 3: (Color online) Transient PM spectra at various time 
delays; arrows indicate the blueshift. The solid lines are fits 
based on Eqs. in the text. Inset: first derivative of Tq 

spectrum; the arrow points to the frequency of the EOT band 
maximum in To spectrum. 



Figure [3] shows several transient PM spectra at vari- 
ous times, t following the pulsed excitation, in the fre- 
quency range that roughly covers the main EOT band. 
Although the data becomes noisier close to the lowest 
Wood's anomaly, a clear deviation from linear behav- 
ior emerges that cannot be reproduced by a simple first 
derivative of the EOT band (Fig. [3] inset); this indicates 
that for plasmonic lattices the various spectral features 
in To{uj) have distinctively different response to a pulsed 
excitation. A noticeable blueshift (~ 70 meV) of the 
main EOT band is apparent that persists for times up to 
100 ps. The blueshift indicates that the EOT is sensitive 
to the imaginary part of permittivity, £^(w), for the fol- 
lowing reason. Within the Drude model for Al, we have 
e'^{uj) = l-Wp/(w^-|-7^). With such and for w < Wp, 
the above-mentioned SPP Bragg condition leads to the 
following expression for, e.g., the (0,1) mode frequency: 



1/2 



(1) 



where tOg = eg is the frequency associated with the re- 
ciprocal lattice period (c is the speed of light). Conse- 
quently, if SPP is the cause of the resonant transmission, 
a transient increase in 7 should lead to a weak redshift of 
the resonance frequency, in disagreement with the data. 
We thus conclude that the observed spectral blueshift 
cannot be explained by the present version of the SPP 
mechanism for the EOT phenomenon [l|. At the same 
time, an increase in 7 causes an increase in the atten- 
uation, e'^, which for 7^0; can be approximated by 
« jujp/uj^. If £^ participates in determining the po- 
sition and lineshape of the EOT bands, then the increase 
in 7 implies higher resonant frequencies (blueshift), in 
agreement with the data. 

For incorporating the non-trivial role of attenuation 
in determining the EOT resonant bands, we adopt an 



analytical effective-medium model [l5| that is a combi- 
nation of the one-mode model 0] and weakly modulated 
permittivity model [l^; these are not directly related to 
SPP excitations. In our model, transmission through the 
film occurs via coupling of the incident light to the eigen- 
modes of the periodic system 8]. Consequently, Tq{u!) at 
normal incidence is given by 



Tn = 



^12^236 



1 + 7'i2r23e ^K/i 



(2) 



where h is the film thickness, ti2 — 2h2/{ni -\- n2) 
and ^23 = 2n3/(n2 -I- ^3) are transmission amplitudes 
through air/metal and metal/glass interfaces, and rki = 
(fik — ni)/{nk + fii) are the corresponding reflection am- 
plitudes, K is the fundamental mode eiegenvalue, and 
fikiuj) are effective refraction indices of each media, which 
are determined by overlap integrals between electromag- 
netic eigenmodes of the perforated film and dielectric 
media The overlap integrals, which are defined at 
the metal/dielectric interfaces, are evaluated using weak 
modulation approximation [iq . whereas the eigenvalue 
K, which is determined by the bulk part of the perfo- 
rated metal, is calculated by using the analogy that ex- 
ists between Maxwell and Schrodinger equations within 
a complex Kronig-Penney potential [l5|. The fit to To{lu) 
using this model with appropriate parameters taken di- 
rectly from the experimental conditions (wp = 15 eV and 
7 = 0.2 eV for Al) is shown in Fig. [TJ The main EOT 
band and first Wood's anomaly in Tq{uj) are accurately 
described by our model; at higher frequencies, where sec- 
ond order features are apparent, the fit is less accurate 
due to the one-mode approximation used in the model. 

To model the ultrafast dynamics at low photoexcited 
carrier density, we calculated the changes in Tq [Eq. ([2])], 
by replacement em(w) em(w) -I- Sem{^,t), where 



dem{uJ,t) = iui, 



■. l[TL{t)] -70 



(3) 



is the photomodulated Drude permittivity. The e-p 
contribution to the relaxation rate in Eq. ^ is given 
by lep{TL) = B{TL/@DfJ[QDlTL), where J{y) = 
/(f dxa;4coth(a;/2)/2, Od « 380 K is the Debye tem- 
perature for Al, and B is a material-dependent constant 
[l7|. Importantly, due to the strong e-p interaction in 
Al, the energy exchange between the hot electrons and 
phonons occurs simultaneously with the thermalization 
of non-equilibrium electron population via e-e scatter- 
ing. This causes the short-time dynamics in Al to be 
an inherently non-equilibrium process. However, since 
the phonons gain energy only from the hot electrons, 
their distribution can be considered as quasi-equilibrium 
[ll,[ll|, and thus characterized by a time-dependent tem- 
perature T^it) that satisfies the rate equation 



Cl-^ = G{t: -Tl) + SoLit), 



(4) 
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where Cl = 900 Jkg ^ is the Al lattice heat capac- 
ity, G — 3.1 X 10^^ Wm^'^K^^ is the e-p interaction con- 
stant, and T* = / dEf{E,t) [1 - fiE,t)] plays the role 
of effective temperature for the athermal non-equilibrium 
electrons having a distribution function f{E, t) (which for 
thermal distribution coincides with the usual hot electron 
temperature) . The last term in Eq. ^ describes the heat 
exchange between the Al film and glass substrate; it is 
given by SGL^t) — kq dTo/dz, where Tq and kq — 1.37 
Jm~^K~"'^ are the glass temperature and thermal con- 
ductivity, respectively. The rate equation ^ is supple- 
mented by heat transport equation in the glass that gov- 
erns the long-time dynamics [l^. In contrast, the short- 
time dynamics is dominated by the energy transfer from 
the non-equilibrium electron population to the lattice, 
and was simulated here by a full numerical solution of 
the Boltzmann equation for f{E,t) [isj . 

The results of our numerical simulations are shown in 
Figs, m and [31 The calculated dynamics are in excellent 
agreement with the data, including the PM rise with sub- 
picosecond time constant [Fig. [H^a)] as well as the long- 
time decay [Fig. (H^b)]. The two insets in Fig. [2] show 
evolutions of the calculated lattice and effective electron 
temperatures. PM dynamics can be summarized as fol- 
lows. After photoexcitation, the effective electron tem- 
perature T* shows a sharp rise that indicates thermal- 
ization of the electron gas via e-e scattering. The sub- 
sequent drop of T* for t < 1 ps after reaching its peak 
value at ~ 200 fs, which is accompanied by an increase 
in Tl and hence in AT/Tq, indicates energy intake by 
the lattice. Note that the energy transfer to the lattice 
takes place while the electron gas is still in the thermal- 
ization stage, in contrast to the dynamics in noble met- 
als [l^ [l^ . This is confirmed by numerical calculation 
of f{E,t) (not shown) indicating the completion of elec- 
tron thermalization by < 1 ps [3. As the electron 
and phonon subsystems equilibrate at tx, PM reaches a 
plateau persisting until the slower process of heat transfer 



to the substrate takes over [Fig. [2Kb)]. 

Fits of the transient PM spectra based on our model 
are shown in Fig. [3] for t > It- The fits capture the 
main EOT band blueshift as well as the nonlinear be- 
havior of the PM spectrum close to the Woods anomaly. 
The nonlinearity stems from the fact that the transmis- 
sion spectrum in plasmonic lattices does not respond uni- 
formly to an increase in the attenuation. This happens 
since the AR frequencies are determined by the lattice 
periodicity, and thus, in contrast to the resonant bands, 
they do not shift when 7 changes. The PM blueshift can 
be traced back to AT dependence on e^, as discussed 
above. The eigenmodes in a plasmonic lattice are the 
result of interference of multiple waves from the perfo- 
rated metal/dielectric interfaces. In the presence of at- 
tenuation, the scattered wave phase has a contribution 
cx Re^/e^ = Em/vprJ ~ l^p/^"^- Therefore, a tran- 
sient change in attenuation causes corresponding shifts 
in the scattering phases. This, in turn, changes the in- 
terference of waves that form an eigenmode. To compen- 
sate for the phase shifts, a transient increase in 7 must 
be accompanied by a corresponding increase in and 
consequently the transmission resonances blueshift. 

In summary, we investigated the ultrafast dynamics 
of nanohole array on Al film that reveals an important 
role played by the plasma attenuation in the makeup 
of surface electromagnetic modes. The transient opti- 
cal spectra discriminate between electromagnetic exci- 
tations with different response to changes in the metal 
dielectric function. The observed transient blueshift fa- 
vors the interference-based mechanisms of explaining the 
EOT phenomenon, as opposed to SPP excitation as the 
cause for the resonant transmission bands. 
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